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Abstract 

A generalized formalism of the so-called non-adiabatic quantum molecular dynamics is presented, 
which applies for atomic many-body systems in external laser fields. The theory treats the nu- 
clear dynamics and electronic transitions simultaneously in a mixed classical-quantum approach. 
Exact, self-consistent equations of motion are derived from the action principle by combining 
time-dependent density functional theory in basis expansion with classical molecular dynamics. 
Structure and properties of the resulting equations of motion as well as the energy and momentum 
balance equations are discussed in detail. Future applications of the formalism are briefly outlined. 
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I. INTRODUCTION 



The non-adiabatic dynamics of electronic and nuclear degrees of freedom in atomic many- 
body systems represents one of the fundamental processes in different areas of physics and 
chemistry. 

Experimentally, exceptional large progress has been made during the last decade in 
studying non-adiabatic processes, in particular in molecules and atomic clusters. So, ex- 
periments with intense femto-second laser pulses interacting with molecules or atomic 
clusters UJ, [| (|, |5], |5], [Tj, |8| have revealed a variety of fascinating new, typical non-adiabatic 
phenomena like the production of keV electrons J2J, MeV ions P| and intense x-rays [Q; the 
Coulomb explosion || connected even with nuclear fusion ||; the multiple plasmon exci- 
tation and relaxation in metallic clusters [0], or the unexpected enhanced ionization with 
decreasing laser intensity ||. Moreover, pump-probe experiments allow now to investigate 
the time-resolved non-adiabatic dynamics, e.g. of photoinduced isomerization processes (for 
a review see @). Finally, refined scattering experiments involving metal clusters |10| and 



fullerenes |TT| revealed detailed insight into electronic and vibronic excitation mechanisms, 
as well as their coupling and related fragmentation processes in those complex systems. 

Theoretically, the non-adiabatic coupling of electronic and nuclear dynamics is one of the 
most challenging problems of atomic many-body theory and, in principle, requires the solu- 
tion of the full time-dependent electron- nuclear Schrodinger equation. At present, however, 
a full-scale numerical solution is barely feasible for the smallest possible molecular system, 
the H^~molecule \T2\. Thus, for larger systems like atomic clusters, phenomenological models, 
based on classical mechanics and/or hydrodynamics |13|, |H] , 15| , |16[ [17], |18| have been devel- 
oped to investigate the mechanism of the intense laser-cluster interaction. More microscopic 
approaches are based on electronic time-dependent Thomas-Fermi theory |19|, or 
related semiclassical (meanfield) approximations |^2|, |23| coupled with molecular dynamics 
(MD) for the nuclear motion. The most advanced microscopic theory to study the coupled 
electronic and ionic dynamics in intense laser-cluster interaction developed so far, is based 
on time- dependent (TD) density functional theory (DFT) in local density approximation 
(LDA) for the treatment of the electronic system coupled with classical MD for the nuclear 
(ionic) dynamics [[24], ^5j. In this approach, the TD-Kohn-Sham equations are numerically 



solved on a grid with the consequence that full 3D calculations [23| are still on the edge 
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of available computational facilities. Therefore, the upper most applications of this theory 
have been obtained within an effective two-dimensional approximation (see |2Bj for a 
review) . 

An alternative fully microscopic approach to the nonadiabatic dynamics in atomic many- 
body systems is the so-called nonadiabatic quantum molecular dynamics (NA-QMD), devel- 
oped recently |27| . In this method, electronic and vibrational degrees of freedom are treated 
simultaneously and self-consistently by combining classical MD with TD-DFT in a finite- 
basis-set expansion of the Kohn-Sham-orbitals. The formalism [^7j has been worked out 
for conservative systems, in particular to investigate adiabatic and non-adiabatic collisions 
involving molecules and atomic clusters. So the NA-QMD theory has been successfully 
applied so far for the description and interpretation of fragment correlations in collision- 
induced dissociation [pgfl , charge transfer cross sections [^. |30|. as well as the excitation 



and fragmentation processes |32], [3l| in collisions of atoms (ions) with small sodium clusters 
and systems as large as fullerenes. 

In this work, we present a generalization of the NA-QMD formalism p7| (hereafter refered 
to as I), suitable to describe also the interaction of large, but still finite atomic many-body 
systems with external laser fields. We derive and discuss the exact equations of motion in 
a systematic way, starting from a general mixed classical-quantum treatment. Energy and 
momentum balance equations are derived as well. Necessary approximations and possible 
simplifications to the exact equations of motion as well as future applications of the formalism 
are briefly summarized. 



II. THEORY 



A. General mixed classical-quantum treatment 

We consider first the general case of a mixed classical-quantum system consisting of 
interacting particles. The Aj classical particles are described by their trajectories R = 
{Ri(t), . . . , RjV;(t)}. They are determined by an explicit time-dependent potential U(H, t) 
as well as the interaction with a system of N c quantummechanical particles, described by 
their many-body wave function \l> = \P(ri, . . . , rjv e , t) (We omit the spin index). This is 
determined by an explicit time-dependent Hamiltonian H(R,, t) which on his part depends 
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parametrically on R. The action of such a system consists of a classical and a quantum part 

A = A C + A q (1) 

with 



A r 



and (atomic units H = e = m e 



A„ 



* 1 ( N- ~) 
to { A > 



(2) 
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1 are used) 
d 



dt 



H(R,t) 



V)dt 



(3) 



with M A the mass of the classical particles and the brackets (...) denote integration over 
all coordinates ri, . . . , rjy e . The equations of motion for the trajectories R and the many 
body state are obtained by making the total action stationary, leading to 



6A 



_6A_ 
SR A (t) 








i^\*)=A(R,t)\9) 



M A R A 
A 











dR, 



■H(R, t) 



(4) 
(5) 



dR A 
l,...,Ni . 

Equations (f|) and (^|) have to be solved simultaneously. They represent the general equations 
of motion of the mixed classical-quantum system defined above. They are much more 
universal than those derived in I from energy conservation. Here they are obtained from a 
general action principle where both, the potential U (defining the classical system) and the 
Hamiltonian H (defining the quantum system as well as the coupling to the classical one) 
may explicitely depend on time. There is no energy or momentum conservation, nevertheless 
classical motion R(t) and quantum dynamics are coupled self-consistently owing to 

the action principle. 

In the next subsection, the potential U and the Hamiltonian H will be specified for an 
atomic many body system, we are interested in. 



B. Atomic many body system 

Considering now N\ ions (nuclei) with charge Z A (A = 1, . . . , N{) and N e electrons exposed 
to an external laser potential (usually, but not necessarily, described in dipole approximation 
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Vl(x, t) = —x ■ E(£), with E(£) the electric field strength) the potential energy of the nuclei 
reads 



U ^ *) = E |r/_r h | - E Z aVl(Ka, t) (6) 
and the Hamiltonian becomes 

h(r, t) = J2^ + J2 v ^ R ^) + E r^—\ (7) 



,-i.A — -tV,R 



r,j — r,- 

i=l i=l i<j J 



with the single particle kinetic energy operator i = The external single particle poten- 
tial V(r, R, t) contains the laser field V^(r,t) and the electron-ion interaction Vi nt (r, R) 



V(r, R, t) = V int (r, R) + V L (r, t) (8a) 
= + • (8b) 

A=l 1 1 

The first term in (8) is time dependent via R(t) and the second one explicitely depends on 
time. Using the definition of the single particle density 

p(r, t)=N c - J d 3 r 2 . . . d 3 r Ni ^*(r, r 2 , . . . , r Nc , t) ■ *(r, r 2 , . . . , r^, t) (9) 

it becomes apparent that the general Newton-type equation ([5]) drastically simplifies with 
(0), (8), (§ leading to 

M^ = -^-/Ap( r , t )«) A=l,...,« . (10, 

Thus, the electronic forces acting on the nuclei are determined by the single particle density 
p(r, t) alone, which is the key quantity in DFT. So, in the next subsection we will reformulate 
the whole problem using TD-DFT to describe the electronic system. 



C. TD-DFT coupled with MD 

According to the basic theorems of TD-DFT |33| any observable of a many body system 
can be expressed as functional of the single particle density (P) and this density can be 



5 



obtained from a non-interacting reference system according to the ansatz 



3=1 



(11) 



with \P J '(r, t) the time dependent Kohn Sham functions. The quantum mechanical part of 
the action (]!]) now reads 



-4 



to 



.9 „ 



pot 



(12) 



where the brackets (...) = j v d 3 r denote integration over the single particle coordinate. 
The potential part in flT2]) 



ti 



.4 



pot 




p(r,t) h/(r, R,i) + 



1 /" p(r/ '|cZV)cZ 3 rcZt + A xc [p] 



r — r 



(13) 



/o 



is a functional of the density p(r,t) and contains the exchange-correlation contribution A xc . 
In concrete applications of TD-DFT, the latter is subject of adequate approximations, like 
the time dependent local density approximation (TD-LDA) or the time dependent optimized 
potential method . In this paper we will not specify A xc and, thus, are dealing with exact 
equations of motion. 

In this sense, variation of (|T^), (0) with respect to the KS-orbitals leads to 



5A 



whereas, variation of (Q), (|l|) and (|T3"D with respect to the trajectories gives 



(14) 



In 



5A 



M A R A 







Ac 



3=1 







dR A 
A=l,...,Ni 

the effective single particle potential V e g(r,R,t) is defined as 



M„ t (r,R) 



(15) 



l/ cff (r,K,t) - — — — 
6p[r,t) 



= V(r,R,t) + J 



P( r ', t) ^ 6A xc [p] 
|r — r'| 5p(r,t) 



(16) 
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In (|15|), the interaction potential Vj nt (r, R), as part of V(r, R, t), is defined according to (8). 

The resulting equations of motion (|I4D, (p~5|) are completely equivalent to (f|), (|5]) and 
accordingly to (|10D. So, with the help of ( |Tl| ) one immediately realizes that (|TJ|) is identi- 
cal to The many body Schrodinger equation ([|), however, is now replaced by a set 
of N e coupled integro-differential single particle KS-equations (|TJj). In the present form, 
these equations have to be solved numerically on a grid, which still is very demanding (if 
not impossible, at present, for large systems in intense laser fields; see also discussion in 
the next section.). A drastic simplification can be achieved, if the (3+l)-dimensional KS- 
orbitals ^(r, t) are represented in a finite basis set. This, however, complicates the classical 



equations of motion (|15D considerably as will be discussed in the next section. 



D. TD-DFT in basis expansion coupled with MD 

In this section, we derive the final equations of motion of the general NA-QMD formal- 
ism and discuss their properties, in particular the resulting energy and momentum balance 
equations. 

The central starting point is to expand the time dependent KS-orbitals ^(r, t) in a local 
basis {<p a } 

vFM) = ^<(i)0 Q (r-R A J (17) 

a 

with the expansion coefficients a J a (t) and the symbol A a denotes the atom to which the 
atomic orbital <p a is attached. 

Although technical details are not the topic of this paper, we note in passing, that the 
use of the linear combination of atomic orbitals (LCAO-ansatz ([H])) has clear advantages as 
compared to a direct numerical solution of the Kohn-Sham equations (|14|). First of all (and 
obviously), the (3+l)-dimensional problem(|l4])) will be reduced to a one-dimensional one 
for the determination of the coefficients a J a (t). Second (and especially important, if intense 
laser fields are considered), electrons with basically different spatial extensions (strongly 
bound core electrons, binding valence electrons as well as practically free electrons in the 
continuum) can be naturally included in the dynamical treatment, provided appropriate 
basis functions <fr a are taken into account . 



With the ansatz flTTj ) the explicit expression of the density is given by 



J'=l a/3 



(18) 



Owing to the implicit time-dependence of the basis a (r— R^ a ), the partial time derivative 



in the action (0) has to be replaced by 



d_ 
dt 



d _ d d 



dt dt ^ dm a 

A=l A 



(19) 



For the following considerations it is convenient to introduce the following matrices: 
the kinetic energy matrix 

Tap := (<p a \i\ 4>p) , 

the hamilton matrix 

H a (3 := (4> a \ t + V e s | 
containing the effective potential V e s defined in (|16]), the overlap matrix 



S a (3 '■= (4>a |0/3) 



the non-adiabatic coupling matrices 

B a (3 '■ = 

which due to (Gil) contains the vector matrices 



d_ 

dt 



B a(3 :z 



and finally, the double differential matrix 







dR, 



C A - 

^a(3 



d 
It 







In addition, we define the transposed matrices 

B 



'a(3 ■' 



D a(3 



Vy a(3 



d_ 

~dt 



D 







dR A 

d 
dR~A 



d 
~dt 



/3a 



BA* 
(3a 



CA* 
(3a 



(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
(27) 
(28) 
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With these definitions and the ansatz flT7|) the quantum mechanical action (|T2]) can be 
written as 



11 



A q = J F q (t)dt - A pot 



(29) 



to 



with 



j=l a/3 



(30) 



The final equations of motion are now obtained by independent variation of the total 
action with respect to a J a (t) and R A (t). With 



SA OF, 



I d 3 r d P SA Pot 



Sa J a{t) dai* J ~ ' dai*Sp(r,t) 
this yields the Kohn-Sham equations in basis representation 







(31) 



(32) 



Pi 



and using Euler's equations 

SA 



dF q d dF q SA pot 



6A r 



SR A (t) dR A dtdR A SR A (t) SR A (t) 
one obtains after some algebra the classical equations of motion 



(33) 



,3* 



OR, 



j=l a/3 



dH, 



+ aj A = 1,...,M (34) 



with the matrix 



3 



(Kfi-v) 



n A+ f^A i Wr+ o-lni -r>A+ q-1 r ~\ 
_ + l^'V _ ^"7 *V "°^J 

7<5 



• (35) 



Equations (|3~2"D and ([34]) represent the central result in the derivation of the generalized 
formalism of the NA-QMD. The Kohn-Sham equations fl32|) are formally very similar to that 
derived in I, but now also contain the laser field which couples to the electronic system (see 
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definition H a p (|21f) , V e g (|16l), V(8) ). The classical equations ( |34"D include the laser field as 
well, which here acts on the nuclei (see definition of U(R,t) (|6|) ). Moreover, the quantum 
part of the forces in ( |34"D differs appreciable from that derived in I obtained from energy 
conservation. In particular, the last term in ( P^D results from the variational principle. It 
represents an important contribution if the momentum balance is considered in the present 
basis set formalism (see below). Obviously, this term vanishes if the basis is complete, i.e. if 



a/3 



(36) 



holds. It will be shown below, that in this case also the remaining terms of the electronic 
contribution to the forces in (0), ([35]) are drastically simplified. In any practical applications 
of the formalism, however, the completeness relation fl36|) can never be fulfilled, and thus, 
the full equations of motion (R3j) have to be considered. 

At first glance, the complicated structure of the forces in ([£§), fl3"o| ) makes it difficult to 
give a transparent interpretation of the correction term resulting from the basis. From the 
theoretical point of view it is therefore very useful to present ([34]), (|35|) in an alternative 
(operator) form and rederive the KS-equations Cj32|) from a basis constrained single particle 
hamiltonian defined as 

h' = t + V eS + X (37) 
with h = i + V e g the usual KS-hamiltonian from (|T4"D, ( |T5] ) and the additional operator 

X := PhP-h + i(l- P)B -iB + (l- P) (38) 
defined with the projectors 

(39) 



a(3 



and 



a/3 



d_ 

~d~t 



o-l 



(40) 



Obviously X vanishes for a complete basis (p6|) . 

With (|38|), fl39l), (|40|) the classical equations of motion (34) can now be rewritten as 

Ne I r) f) 

h ) 4 ( 41 ) 



<9R, 



j=l a{3 



d d « 

K nt (r,R) + ^X 



dR 
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leading finally, with (|T7|), to 

dU(R,t) 



M A R A = - 



(42) 



In addition, the equations of motion (^) are equivalent to the standard form of the time- 
dependent KS-equations 



= (t + V eS (r,R,t) + X)V 

however, with the additional single particle operator X (|38|). 
inserting the ansatz ( O ) into ( |43"D which leads to 



(43) 

This can easily be seen by 







(44) 



and, therefore, finally to (^) because the basis {4> a } must be linearly independent. 

The implicit equations of motion fl4"2"| ) and p3] ) are thus completely equivalent to the 

), used in practical calculations. They allow however a more 



explicit expressions (|32|) , 
transparent interpretation of the present theory: 

The use of a finite basis expansion has the same effect as the introduction of an additional 
operator in the hamiltonian. This is similar to the introduction of constraining forces in 
classical mechanics, if the dynamics is investigated under boundary conditions. Further, one 
can now explicitely see that the "coupled channel" equations (^2]) and the "constrained" 
forces (|34]), (|35|) reduce to the standard KS-equations ( |i~4"D and Newton-equations (|i~5D, 
respectively (|T0|), if the basis is complete. 



E. Energy and momentum balance 



In order to derive the energy balance we define an exchange-correlation energy according 



to 



with the important property 



A xc [p] = J E xc [p](t)dt 
to 



SA xc [p) = 6E xc [p](t) 
6p(r,t) 6p(r) 



(45) 



(46) 
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Note the different arguments in dp on the left and right hand side of fl4q). With this, the 
potential energy of the quantum system can be written as 

E pot (t) = J p(r,t)(v(r,R,t)+ 1 - j t^^A d 3 r + £ xc [p](t). (47) 



d_ 
~dt 



The total time derivative of this functional is given by 

^ P ot(t) = J d 3 rj t p(r,t)-V cfi (r,K,t) + j d 3 r p(r,t) ■ jV(r,R,t) (48) 
with 

d d Ni d 

p M )= p(r,«)+£R ,(,,,) (49) 

A=l A 

and 



|y (r , H , = ^ + £ > ^, W r,R) . (50) 

A=l ^ 



Now, the total energy of the system can be defined 



E(t) = J2^A + u(n,t) + Y,J2< T ^ + E ^[p](t) ■ (si) 

A=l j=l a/3 



The total time derivative of this quantity is obtained after a longer calculation using (p8|). 
(^9|), (p0|) and the equations of motion (^) as 

J A=l 

As expected, this quantity is conserved for vanishing or time-independent external fields. 

A more transparent expression for the energy balance can be obtained in dipole approx- 
imation (i.e. Vl(x, t) = — x ■ E(t) ) leading to 

^=-de(t)E(t) + di(t)E(t) (53) 

with the dipole moments of the electrons 

d e (t) = J p(v,t)vd 3 r (54) 



and the ions 



d i (t) = J^Z A R A (t) . (55) 

A=l 
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From this expression it is clearly seen, that in a homonuclear system (Za = Z =const.) the 
ions will not be excited by the laser, because in the center of mass system the nuclear dipole 
moment vanishes, i.e. 



di = Z r a = . (56) 

A=l 

In order to obtain more insight into the electronic excitation (deexcitation) process, it is 
convenient to consider the total energy change 



AE A = - d c (t)E(t)dt (57) 

J —oo 

together with the Fourier-transformed dipole moment 

d c (u) = — / e-^d e (t)dt . (58) 

One now immediately realizes that in a continuous wave field 

E = 3?(E e~ iaJL *) (59) 

the electronic system adsorbs (desorbs) energy only if the imaginary part of d c (u) does not 
vanish at the laser frequency to = u;l, i.e. 

AE cl = u L %(d c (uj L ) -E ) . (60) 

In the linear response region this is the case only if uj^ coincides with the excitation energy 
of an optical excited state. 

In the other extreme case of very short laser pulses 

E = E 5(t) (61) 

all frequencies do contribute simultaneously to the excitation (deexcitation) process, i.e. 

AE cl = / dwwS(dc(w) • E ) . (62) 



For finite laser pulses, the total electronic energy change can be obtained by solving (| 
numerically, together with the the full equations of motion (|32"1), fl34|) to calculate the dipole 
moment d c (t). 

We note also, that the present formalism can be favourably used to calculate optically 
excited states (i.e. Born-Oppenheimer surfaces) as well as optical excitation spectra in the 
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linear response region from ( |58j) by solving the KS-equations ( p2[) for fixed nuclear position R 
and "numerically short" 5-pulses (^1|). Details of this procedure will be discussed elsewhere 

In order to investigate the momentum balance we start with the total momentum 

P = P c + P q (63) 



as the sum of the classical part 



P C = J^M A R A (64) 

A=l 



and the quantummechanical part 

N a Nr. 



j=l j=l a/3 



Using the identity 

. a (r - K Aa ) = -g^-M* ~ H*0 = - E MT^ r ~ ( 66 ) 

A a A=l A 

the latter one can be transformed into 

j=l a/3 A 

Now, the total derivative with respect to time can be obtained using fl3"2|) and ([34]) leading 
after an extensive calculation to 

-P = - J p(r, t) V(y L (r, t) + Kc(r, t))rf 3 r + ^ Z A W l (Ra, *) • (68) 

Besides the expected dependence on the laser field, this balance contains a term that depends 
on the exchange-correlation potential V xc = 5A XC /Sp(r, t). This one vanishes for the exact 
Vxc, which is a general property of TD-DFT ||38|| . Without this term one also immediately 
realizes, that in dipole approximation the total momentum balance vanishes for neutral 
systems, i.e. 

TV; 



dt 



I p(r, t)d 3 r + Y J Z A \ E(t) = (69) 

J A=l J 
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which is due to the classical, not quantized treatment of the laser field. 

We note finally that the momentum balance (|68| ) can be derived also (and much simpler) 
without basis expansion. The derivation, carried out here, therefore proofs nicely the validity 
and stresses the importance of the finite basis correction terms in the forces (|53fl, (PSj) 
following from the variational principle. 



III. SUMMARY AND OUTLOOK 



We have derived in a systematic way a generalized formalism of the NA-QMD which 
applies for finite atomic many-body systems in external fields. It is based on a mixed 
classical-quantum approach where the electronic system is described by TD-DFT in local 
basis expansion and the nuclear degrees of freedom are treated classically by molecular 
dynamics. Self-consistent equations of motion are derived from a general action principle. 

We have presented here the exact equations of motion. They can be solved without 
further approximations for one electron systems, like or HD + []36| where the exact 
exchange-correlation term cancels the Hartree-term in the effective potential flTEp. For many- 
electron systems, approximate equations of motion, as derived e.g. in I on a tight-binding 
level, can be obtained from the general formalism as well. We intend however, to realize 
the numerical implementation of the whole formalism also on the ab-initio level using the 



time-dependent optimized potential method [35| for the exchange correlation part in the 
action ([13]). Preliminary results obtained within this method for organic molecules, like 
ethylene C2H4, show excellent agreement with Cl-calculations E9] concerning the ground- 



state properties (i.e. bonding lengths, angles etc.) as well as optical excitation spectra ||0 



As a first application of the whole time-dependent formalism we intend to investigate the 
cis-trans isomerization process of C2H4 in short laser pulses pEOfl . 

Another very interesting and fascinating field of application concerns the excita- 
tion, ionization and fragmentation mechanism of atomic clusters in intense laser fields 
0, |3|, [|, |^, ^ 0, [|. Here an all electron treatment together with an appropriate descrip- 
tion of the continuum in the ansatz ( |17|) is required which, as discussed in the text, can be 
incorporated in the present formalism p7| . 

This work was supported by the DFG through Forschergruppe "Nanostrukturierte Funk- 
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